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ABSTRACT 

Current X-ray observatories make it possible to follow the evolution of transient and variable X-ray binaries 
across a broad range in luminosity and source behavior. In such studies, it can be unclear whether evolution in 
the low energy portion of the spectrum should be attributed to evolution in the source, or instead to evolution 
in neutral photoelectric absorption. Dispersive spectrometers make it possible to address this problem. We have 
analyzed a small but diverse set of X-ray binaries observed with the Chandra High Energy Transmission Grating 
Spectrometer across a range in luminosity and different spectral states. The column density in individual photo- 
electric absorption edges remains constant with luminosity, both within and across source spectral states. This 
finding suggests that absorption in the interstellar medium strongly dominates the neutral column density ob- 
served in spectra of X-ray binaries. Consequently, evolution in the low energy spectrum of X-ray binaries should 
properly be attributed to evolution in the source spectrum. We discuss our results in the context of X-ray binary 
spectroscopy with current and future X-ray missions. 



constant or jointly determined, this method is adopted some- 
what arbitrarily, without specific reference to an observational 
or theoretical basis (see, e.g., Sobczak et al. 2000; Miller et al. 
2001; Nowak, Wilms, & Dove 2002; Park et al. 2004; Gierlin- 
ski & Done 2004; Kalemci et al. 2005; Rykoff et al. 2007). 

Dispersive X-ray spectroscopy can separate variations in ab- 
sorption from variations in the source spectrum. Even the reso- 
lution afforded by CCD spectrometers blends many weak edges 
into the continuum, and the optical depth and wavelength of 
stronger edges can be difficult to measure precisely. For in- 
stance, in fits to CCD spectra, the multiple Fe L edges can be 
modeled satisfactorily in terms of a single step function. In con- 
strast, Chandra High Energy Transmission Grating Spectrome- 
ter observations of X-ray binaries have resolved the individual 
Fe L edges and their detailed properties (see, e.g., Schulz et al. 
2002, Juett et al. 2006). Beyond simple edges, dispersive X-ray 
spectroscopy has clearly revealed the neutral O I ls-2p reso- 
nance absorption line that was first glimpsed with the crystal 
spectrometer aboard Einstein (Schattenberg & Canizares 1986, 
Paerels et al. 2001). 

In this Letter, we extend recent work using high resolution 
X-ray spectroscopy by examining the evolution of neutral ab- 
sorption in a sample of X-ray binaries. We find no evidence 
for variable absorption, consistent with an interstellar origin for 
neutral absorption. 

2. DATA SELECTION AND REDUCTION 

Our analysis is restricted to observations made with the 
Chandra High Energy Transmission Grating Spectrometer. The 
HETGS offers better resolution than the XMM-Newton Reflec- 
tion Grating Spectrometer, as well as coverage across the full 
0.3-10.0 keV band. When the equivalent neutral hyrdogen den- 
sity is too high, much of the source continuum spectrum is 
scattered out of the line of sight. This must be avoided as 
absorption features are measured relative to a continuum. A 
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1. INTRODUCTION 

X-ray binaries show a remarkable variety of extreme phe- 
nomena, including rapid variability in many wavelength bands, 
the production of relativistic jets and hot winds, and outbursts 
that span factors of 10 6 in luminosity (for a review of black hole 
transients, see Remillard & McClintock 2006). Monitoring ob- 
servations with the Rossi X-ray Timing Explorer and Swift en- 
able the evolution of X-ray spectral and timing properties of 
X-ray binaries to be traced in order to understand the accretion 
process (see, e.g., Rykoff et al. 2007). Model degeneracies that 
may occur in fitting a single spectrum can often be resolved by 
examining how multiple spectra evolve with time. 

Major changes within the outbursts of transient and persistent 
X-ray binaries are typically linked with changes in the nature 
of the accretion disk (see, e.g., Esin, McClintock, & Narayan 
1997; Homan et al. 2001; also see Miller et al. 2006a). Track- 
ing accretion disk evolution is therefore a primary goal of mon- 
itoring observations. The disk is directly visible in soft X-rays 
in black hole X-ray binaries, and likely also in neutron star bi- 
naries (e.g. Cackett et al. 2008a). However, the K shell pho- 
toelectric absorption edges from most abundant elements fall 
below 2 keV, as do Fe L shell edges. This absorption can com- 
plicate efforts to study disk properties, especially when using a 
spectrometer with limited energy range and resolution. 

Prominent theoretical treatments of accretion flows predict 
modest winds (e.g. Begelman, McKee, & Shields 1983) and/or 
outflows with a high ionization parameter (e.g. Narayan & Ray- 
mond 1999). Observationally, even the most extreme winds de- 
tected in X-ray binaries do not appear to generate enough col- 
umn density to add appreciably to the neutral absorption col- 
umn (Miller et al. 2006). Nevertheless, it is not uncommon to 
allow the equivalent neutral hydrogen column density to vary in 
modeling multiple spectra of an X-ray binary (see, e.g. Brock- 
sopp et al. 2006; Caballero-Garcia et al. 2009; Cabanac et al. 
2009). Even in cases where the neutral column density is held 
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number of clear absorption edges are detected in Cygnus X- 
1, which has an equivalent neutral hydrogen column density 
of N H = 6.2 x 10 21 cm -2 (Schulz et al. 2002). However, an 
estimated column density of 8.0 x 10 21 cm" 2 does not allow 
the absorption edges in XTE J 1550-564 to be studied in de- 
tail (Miller et al. 2003). We therefore restricted our analysis to 
sources equal to or below approximately 6 x 10 21 cm -2 . 

Our analysis was also restricted by the need to select vari- 
able and transient sources. Observations of transient sources 
with large missions can be difficult, and relatively few sources 
with low or moderate column densities have been observed on 
multiple occasions. The sources and observations selected for 
our analysis are listed in Table 1 . Different source types, com- 
panion types, and orbital periods are included in the sample. 

When observing bright sources using the HETGS, the ACIS- 
S array of CCDs must be run in "continuous clocking" mode 
to prevent photon pile-up. For consistency, the observations 
considered in this work employed this observing mode. The 
only exception is a single observation of GX 339-4 in a low 
flux state (obsid 4420; see Table 1). An observation of Cygnus 
X-l in the low/hard state that was not taken in "continuous 
clocking" mode (obsid 107) suffered heavy photon pile-up, 
complicating absorption spectroscopy (see Juett, Schulz, & 
Chakrabarty 2004), and was therefore rejected. For more in- 
formation on the specific nature of HETGS observations made 
in "continuous clocking" mode, please see Miller et al. (2006b) 
or Cackett et al. (2008b). 

Most of the spectra considered in this work were obtained 
using the "TGCat" data center (see http://tgcat.mit.edu). This 
facility is run by the Chandra X-ray Center and provides first- 
order HETGS spectra and responses derived using up-to-date 
CI AO software and calibrations. An error in the instrument 
offsets used to observe XTE J1817-330 required custom pro- 
cessing as per the procedure outlined in Miller et al. (2006b) 
and Cackett et al. (2008b). Two observations of 4U 1820-30 
were previously reduced and analyzed by Cackett et al. (2008b) 
(obsid 6633 and 6634), and the same spectra are used in this 
work. Using CIAO version 4.0.2, we combined the first-order 
MEG spectra and instrument response files using the CIAO tool 
"add_gratings_spectra". The same procedure was repeated for 
the first-order HEG spectra. Again owing to the pointing error, 
only the minus-side first-order spectra of XTE J1817-330 were 
used in our analysis. To create distinct spectra of Cygnus X-2 
in the flaring, horizontal, and normal branches of its "Z" track, 
we employed exactly the same selection criteria and procedure 
as Schulz et al. (2009) to both observations (obtaining a total 
of six spectra). The "aglc" script was used in the ISIS analy- 
sis package (Houck & Denicola 2000) to create lightcurves and 
color-color diagrams that were used to create spectra from each 
branch. 

All of the spectral fits reported in this work were made using 
XSPEC version 12 (Arnaud & Dorman 2000). Unless other- 
wise stated, errors are la errors. 



3. ANALYSIS AND RESULTS 

The true neutral absorption column density observed in a 
given spectrum should not be a function of the continuum flux 
model. However, to ensure model-independent results, we fit 
for prominent neutral photoelectric absorption edges 
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FlG. 1. — In the plot above, high/soft (black) and low/hard (red) spectra of 
Cygnus X- 1 in the region of the Ne K edge are plotted as a ratio to a power-law 
continuum. The best-fit column density for the edge has been set to zero to 
show the remarkable consistency of the edge in different spectral states. 

individually - in narrow wavelength ranges - using local power- 
law models for the continuum. The O K, Fe L, and Ne K edges 
were fit in the 20-25A, 16-18A, and 13-15A bands, respec- 
tively. At longer wavelengths, the MEG has a higher effective 
area than the HEG; moreover, it covers longer wavelengths than 
the HEG. Fits to absorption edges were therefore made to the 
combined first-order MEG spectrum from each observation. 

In many of the spectra, it was quickly apparent that simple 
step-function edge models do not fully describe the data. We 
therefore used the "tbnew" model (Wilms et al. 2009), which 
includes up-to-date cross-sections and detailed edge structure. 
A new functionality in "tbnew" can be used to fit the column 
density in an individual edge by setting the overall column den- 
sity and abundance of all elements (apart from the element of 
interest) to zero, and restricting the column density to negative 
values for the parameter of interest. (Negative values instruct 
the model to return the column density for the element; positive 
values return the abundance of a given element relative to hy- 
drogen.) All fitting parameters relating to grain properties were 
fixed at their default value. In our fits, then, "tbnew" functioned 
as a model with only one variable parameter (the column den- 
sity of the element of interest). Due to a combination of a rela- 
tively high total absorption column and the particular choice of 
instrumental pointing offsets used, useful constraints on the O 
K edge could not be obtained from existing data on GX 339-4, 
4U 1820-30, and Cygnus X-2. We note that Yao et al. (2009) 
were able to study the Oxygen edge in Cygnus X-2 by adding 
multiple spectra, but this procedure is inconsistent with the aim 
of this study. 

To estimate the luminosity of each source in a given obser- 
vation, we fit the first-order HEG spectra in the 1.2-10.0 keV 
range using simple absorbed continuum models. The equiva- 
lent neutral hydrogen column density was allowed to vary but 
values were found to be broadly consistent with fits to individ- 
ual edges. These models should be regarded as fiducial. Some 
of the fits obtained are not formally acceptable, owing to cali- 
bration uncertainties and/or unrelated ionized absorption lines 
in the spectra. The black hole sources were fit with a continuum 
consisting of a disk blackbody and a power-law. The neutron 
star spectra were all fit with a continuum consisting of a simple 
blackbody function and a power-law. To represent the evolution 
of absorption with luminosity, the luminosity in each observa- 
tion was calculated for a given distance. Errors on the luminos- 



Miller, Cackett, & Reis 



3 



ity were then calculated using the la errors on the continuum 
flux. The distances assumed in order to calculate luminosity 
values are listed in Table 1 . 

Some of the spectra modeled in this paper have not been 
analyzed previously. All four Chandra observations of XTE 
J1817-330 were made in the high/soft state. This is evi- 
dent from inner disk color temperatures that range between 
kT=0.48 keV and 0.94 keV, and power-law indices that range 
between T = 2.4 and 4.0. The evolution of the bright phase 
of the outburst of XTE J1817-330 is detailed in Rykoff et al. 
(2007), and comparisons to that work again indicates that the 
Chandra observations were all made in high/soft states. Simi- 
larly, spectral results have not been reported based on fits to the 
second and third observations of GX 339-4 that are considered 
in this work (obsids 4569 and 4570). Here again, the flux was 
strongly dominated by the accretion disk, with inner disk color 
temperatures of kT = 0.68 keV and 0.81 keV. Finally, the third 
observation of 4U 1 820-30 (obsid 7032) has not been published 
previously, but analysis as per Cackett et al. (2008b, 2009) sug- 
gests that the source was observed in the "banana branch". 

We find that neutral absorption is remarkably constant with 
source luminosity, and consistent across different spectral 
states. Figure 1 shows the Ne K edge in the high/soft and 
low/hard state spectra of Cygnus X-l. The jump in the contin- 
uum at the Ne K edge is consistent. Figure 2 plots the column 
density measured in individual absorption edges versus lumi- 
nosity for the sources in our sample. For the black hole sources, 
the column density in a given edge is constant with luminosity 
within la errors. The case is much the same with the neutron 
star X-ray binaries. The column density measured in the Ne K 
edge of 4U 1820-30 is marginally different at the la level in 
two observations, but both are consistent with a third observa- 
tion. All three columns agree at the 90% confidence level. The 
same holds true for two observations of Cygnus X-2. 

The Ne column densities that we have measured agree with 
values reported by Juett, Schulz, & Chakrabarty (2004) for 4U 
1820-30, Cygnus X-2, and GX 339-4. The Ne column density 
that we measured in Cygnus X-2 also agrees with a detailed 
analysis made by Yao et al. (2009). Our value of the Ne column 
density agrees with that measured in Cygnus X-l by Hanke et 
al. (2009) though our values of the O and Fe columns differ by 
factor of approximately two. These disparities can likely the 
result of different continuum modeling procedures and imple- 
mentations of "tbnew". 

4. DISCUSSION AND CONCLUSIONS 

To better understand the evolution of the low energy spec- 
trum of X-ray binaries, we made fits to individual photoelectric 
absorption edges in high resolution X-ray spectra of selected 
sources. The column density measured in individual edges is 
not observed to vary across different spectral states, nor over 
a broad range in luminosity (see Table 1 and Figure 2). This 
suggests that gas from X-ray binaries is not typically an impor- 
tant source of the neutral absorption observed in the spectra of 
these systems. Rather, neutral absorption must be dominated 
by the ISM. A similar conclusion was reached by Juett, Schulz, 
& Chakrabarty (2004) based on upper limits on the velocity 
dispersion of the ISM as measured through absorption lines. 
Evolution in the low energy spectrum of typical X-ray binaries, 
then, is best attributed to evolution in the source continuum. 

Neutral absorption in X-ray spectra is often fit by a single 



model that parameterizes the accumulated absorption from in- 
dividual edges as an equivalent neutral hydrogen column den- 
sity. Values obtained from high resolution spectra are likely to 
give the best measure of a true equivalent total column den- 
sity. In practice, instrumental problems such as internal scatter- 
ing, carbon build-up from optical blocking filters, and gain drift 
could prevent the adoption of a gratings-derived value for the 
neutral column. In such cases, our results suggest that a single 
value of the equivalent neutral hydrogen column density should 
be used to fit multiple spectra from monitoring observations of 
a given source with a given detector. 

The outflows that are observed in X-ray binaries are highly 
ionized - dominated by He-like and H-like charge states (see, 
e.g., Lee et al. 2002, Miller et al. 2004, Miller et al. 2006b, 
Schulz et al. 2008). An especially dense wind was observed 
in GRO J 1655-40, and even in that case the ionized columns 
observed are insufficient to create strong absorption edges that 
could be mistaken for additional neutral absorption (Miller et 
al. 2006c, 2008). Moreover, in sources such as H 1743-322, 
GRO J1655-40, and GRS 1915+105, a paradigm is emerging 
wherein ionized winds are active in soft, disk-dominated states, 
but absent in spectrally hard states like those that typically hold 
when sources accrete at a low fraction of their Eddington limit 
(Miller et al. 2006b, 2006c, Miller et al. 2008, Neilsen & Lee 
2009). In Figure 2, it is clear that any variation in ionized winds 
across states does not affect the column density measured in 
neutral edges. 

Our results are based on spectra which only reach down 
to approximately 0.01 L E( jd.. Sensitive high-resolution X-ray 
spectra that would permit strong constraints on absorption vari- 
ability have not yet been obtained from sources at lower lu- 
minosity. However, theoretical arguments again point to ion- 
ized outflows that would contribute little to a neutral column. 
Winds from advection-dominated accretion flows are expected 
to be extremely hot (since advective flows are very hot), and 
line spectra should be dominated by He-like and H-like charge 
states (e.g. Narayan & Raymond 1999). Recent observations of 
the stellar-mass black hole V404 Cyg, which accretes at about 
10" 5 L E( jd., are able to rule-out the winds predicted by some 
advective models (Bradley et al. 2007). This further suggests 
that ionized outflows are not likely to contribute significantly to 
line-of-sight absorption, even as sources approach quiescence. 

There are at least two classes of sources where our results 
may not hold in all circumstances. Neutron star X-ray binaries 
known as "dippers" are viewed at high inclinations, and mate- 
rial in the outer disk can block emission from the central engine 
(see, e.g., Diaz Trigo et al. 2006). Within such flux dips, the 
observed neutral absorption column may vary due to the chang- 
ing geometry within the binary. Similarly, winds from massive 
stars are known to be clumpy and to sometimes cause flux dips; 
some of these dips may also cause variations in the observed 
equivalent neutral hydrogen column density (e.g. Balucinska- 
Church et al. 2000). 

We thank Joern Wilms for creating the "tbnew" model used 
in this work, and for guidance on its implementation. We ac- 
knowledge helpful comments from the anonymous referee that 
improved this paper. J.M.M. gratefully acknowledges funding 
from the Chandra Guest Observer program. EMC gratefully 
acknowledges support provided by NASA through the Chandra 
Fellowship Program, grant number PF8-90052. RCR acknowl- 
edges STFC for financial support. 
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Source Name 


P orb (days) 


distance (kpc) 


Companion 


ID 


Start (MJD) 


Exposure (ksec) 


State 


GX 339-4 


1.76" 


8.5" 


low-mass" 


4420 


52715.8 


76.2 


hard/intm.* 










4569 


53239.5 


50.1 


high/soft' 










4570 


53282. 


50.2 


high/soft' 


XTEJ1817-330 


short^ 


8.5" 


low-mass^ 


6615 


53779.1 


50.0 


high/soft' 










6616 


53790.2 


49.6 


high/soft' 










6617 


53809.3 


46.7 


high/soft' 










6618 


53877.6 


51.0 


high/soft' 


Cygnus X- 1 


5.6 C 


2.5* 


09.7 Iab E 


3724 


52485.7 


26.4 


high/soft'" 










3815 


52702.7 


55.9 


low/hard" 


4U 1820-30 


0.007* 


1.6" 


low-mass s 


6633 


53959.9 


25.2 


BB'' 










6634 


54028.2 


25.1 


BB' ; 










7032 


54044.6 


46.3 


BB' 


Cygnus X-2 


9.843' 


1.2 1 


A5-F2' 


8170 


54337.7 


78.7 


NB, HB, FB« 








8599 


54335.2 


70.6 


NB, HB, FB« 



NOTE. — The table above lists the name, orbital period, and the companion type (where known) for the X-ray binaries included in our sample. The Chandra 
observation identification number, observation start time, and exposure time is also given for each observation. Finally, the "state" in which the source was observed 
is provided. The first three sources are black hole candidates, the last two sources are low-mass X-ray binaries harboring neutron stars. 

" Hynes et al. (2003). * Zdziarski et al. (2004). c Torres et al. (2006). d Based on its coordinated near to the Galactic center. e Gies & Bolton (1982). f Bregman et al. 
(1983). « Stella, Priedhorsky, & White (1987). * Kuulkers et al. (2003). ' Cowley, Crampton, & Hutchings (1979). ' Kuulkers et al. (1999). * Miller et al. (2004). ' 
This work. m Feng, Tennant, & Zhang (2003). " Miller et al. (2006a). p "Banana" branch; Cackett et al. (2009). q Normal branch, horizontal branch, flaring branch; 
Schulz et al. (2009). 
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FlG. 2. — The plot above shows the column density measured in specific low energy absorption edges versus luminosity, for black hole and neutron star X-ray 
binaries. Errors on column density are Icr errors. 



